Acoustic phonons and spin relaxation in graphene nanoribbons 
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Phonons are responsible for limiting both the electron mobility and the spin relaxation time in 
solids and provide a mechanism for thermal transport. In view of a possible transistor function as 
well as spintronics applications in graphene nanoribbons, we present a theoretical study of acoustic 
phonons in these nanostructures. Using a two-dimensional continuum model which takes into ac- 
count the monatomic thickness of graphene, we derive Hermitian wave equations and infer phonon 
creation and annihilation operators. We elaborate on two types of boundary configuration, which 
we believe can be realized in experiment: (i) fixed and (ii) free boundaries. The former leads to a 
gapped phonon dispersion relation, which is beneficial for high electron mobilites and long spin life- 
times. The latter exhibits an ungapped dispersion and a finite sound velocity of out-of-plane modes 
at the center of the Brillouin zone. In the limit of negligible boundary effects, bulk-like behavior 
is restored. We also discuss the deformation potential, which in some cases gives the dominant 
contribution to the spin relaxation rate T^^ . 

PACS numbers: 63.22.Rc, 72.20.Dp, 76.60.Es, 81.07.Oj 



I. INTRODUCTION 

Its interesting electronic, mechanical, and thermal 
properties have made graphene a promising candidate 
for a wide range of apphcations, including ballistic tran- 
sistors as well as spintronics and nanoelectromechani- 
cal devices and heat management.^ ^ There are, how- 
ever, a number of challenges: (i) for epitaxial graphene, 
which is desirable for a controlled, large-scale production, 
the strong coupling to a substrate compromises these 
properties, (ii) graphene has no band gap, a handicap 
for typical semiconductor applications, and (iii) acoustic 
phonons limit the carrier mobility relevant for transistor 
functions.* 

The first issue can be overcome by removing substrate 
material from underneath the carbon layer such that a 
trench is formed and the electronic properties of free- 
standing graphene are restored. ^^ ^^ The second chal- 
lenge is met by graphene nanoribbons (GNRs), graphene 
strips with a width at the nanometer scale (e.g., L ^ 
1 ^m, W ~ 30 nm) which can exhibit a band gap.^® 
Combining these advantages, the free-standing GNR ob- 
tained from epitaxial graphene on a trenched substrate 
is a very interesting design that deserves a detailed dis- 
cussion of its phonons. 

In this paper, we use a continuum model to study 
the acoustic phonon properties and displacement fields 
u{r) = {ux, Uy,Uz) (out-of-plane modes shown in Fig. 1) 
of two different types of GNR that we think can be re- 
alized in experiment: (i) extended graphene that covers 
a thin trench, resulting in a GNR parallel to the trench 
and with fixed lateral boundaries. Fig. 2 (a); (ii) a strip 
of graphene that stretches over a wide trench, leading to 
a GNR perpendicular to the trench and with free lateral 
boundaries. Fig. 2 (b). For both setups, we derive the 
low-energy acoustic phonon spectra from a continuum 
model that respects the monatomic structure of graphene 
and write down the quantum mechanical form of these 




FIG. 1: (Color online) Displacement field of out-of-plane 
modes. The dimensionless coordinates are x — x/W, y = 
y/W, and z = z/W. (a),(b) Fixed boundaries. (c),(d) Free 
boundaries. (a),(c) Fundamental mode. (b),(d) First over- 
tone. 

phonons. 

Our results can be probed experimentally via estab- 
lished techniques like electron energy loss spectroscopy or 
Brillouin light scattering. In addition to the electron 
mobility, the phononic behavior is essential for carbon- 
based nanoelectromechanical systems.^ A recent ex- 
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FIG. 2: (Color online) Two nanoribbon configurations where 
graphene (blue) is spanned over a trenched substrate (orange): 
(a) fixed and (b) free boundaries. The coordinate system 
is chosen in such a way that the undeformed ribbon lies in 
the xy plane and that the y and ribbon axes coincide. We 
assume the ribbon length to be much larger than the ribbon 
width (L ^ and parallel lateral boundaries (dark blue) 
at a; = ±VK/2. 



ample where the electron-phonon coupling has been ob- 
served experimentally is the Franck-Condon blockade in 
suspended carbon nanotube quantum dots.^'^ Phonons 
also give rise to spin relaxation within a time Ti, which 
is important for spintronics devices'* 2* "W\e spin-orbit 
interaction admixes different spin states (f, i) and elec- 
tron orbits fc [see Eq. (18) below]. As a consequence, 
the electron-phonon coupling -/Jepc can mediate the Zee- 
man energy g^BB = fiuj, where g is the electron gr-f actor 
in graphene and /is denotes Bohr's magneton, to the 
phonon bath with phonon numbers n^j. The rate for 
spin relaxation via emission of a phonon with energy huj 
is given by 



EPC 



fcf, ?^t^)Ppstatos(?it^) , (1) 



with an explicit dependence on the phonon density of 
states PstatGs- Several mechanisms contribute to T^^ and 
in some cases the deformation potential^^ gives the 
dominant contribution. If the Zeeman energy lies within 
the energy gap of GNR phonons with fixed boundaries 
and if the temperature is sufficiently low, the spin life- 
time obtained from Eq. ( 1 ) diverges due to a vanishing 
density of states. 



II. CONTINUUM MODEL IN 2D 

Low-energy acoustic phonons at the center of the Bril- 
louin zone have a wavelength much larger than atomic 
distances and thus can be derived from continuum me- 
chanics. The carbon atoms in graphene lie within a two- 
dimensional surface and this property is conserved upon 
deformations, making graphene a quasi-two-dimensional 
material in three-dimensional (3D) real space. Conse- 
quently, all components of the displacement field u{r) 



can be nonzero but the components Uiz of the strain ten- 
sor Uik = (diUk + dkUi)/2 vanish identically. While Uxz 
and known to vanish for thin plates in the xy 

plane in general, the monatomic thickness of graphene 
implies that u^^ must vanish as well. With Uiz = 0, the 
elastic Lagrangian density of monolayer graphene is given 
by 



C = T-V 



2 



A 



(2) 



where A = 9^ 
u,n,+Uzz and ' 



dy , the sum convention with ua — 



yy-TLizz aiiu LL^f. ~ u-2.3, -r Li^y -r • • • has bccu uscd, p is the 
surface mass density, and k is the bending rigidity. ^^ '^'^ 
Note that the 3D bulk elastic constants have been re- 
placed by their 2D analogs A = 2/i/Lt3DA3D/(2/i3D + Asd) 
and /i = /i/^3D where h is the plate thickness. The bulk 
and shear moduli are then given as B = A -I- and /z, 
respectively. 

Application of the Euler-Lagrange formalism to the 
functional ( 2 ) leads to the coupled set of differential equa- 
tions for in-plane modes 



/9 iia; = (i? + \i)dl.Ux + [I dyUx + B dxdyUy , 

Pily = {B + p)dyUy + fl dlUy + B dxdyUx ; 



(3) 



which are decoupled from the differential equation for the 
out-of-plane modes, 



pUz 



dlfuz 



(4) 



Assuming nanoribbon alignment with the y axis, fixed 
boundaries are described by 



dxUz 



(in plane), 
(out of plane) 



(5) 
(6) 



dX X = ±W/2; see Fig. 2 (a). While these boundary con- 
ditions hold for both 2D and 3D lattices, we emphasize 
that lattice dimensionality does affect free boundaries. 
For free edges in 2D it is required that, at x = ±.W/2, 



dxUx + adyUy = 

dxUy + dyUx = 

dluz + (2 - a)dxd1uz - 
dluz + crdyUz = 



(in plane), (7) 
(out of plane), (8) 



where the quantity a denotes Poisson's ratio, Fig. 2 (b). 
Together with Young's modulus £ = h'E^D, a relates to 
the bulk and shear moduli as 



B 



2(1 -a) 



2(1 + a) 



(9) 



III. CLASSICAL SOLUTION 

Typically, the length of a graphene nanoribbon exceeds 
its width many times, '^^ "^^ L ^ W, thus allowing for 





FIG. 3: (Color online) Dispersion relations obtained from the 
procedure described in Sec. Ill, The wavenumber q is given 
by g = qW and the frequency uj of in-plane (out-of-plane) 
phonons by Coxy = ujy^ p/'EW (uiz = p/nW^). (a) In- 
plane modes with fixed boundaries, (b) In-plane modes with 
free boundaries, (c) Out-of-plane modes with fixed bound- 
aries, (d) Out-of-plane modes with free boundaries, (e)- 
(h) Dispersion relations (a)-(d) at the center of the Brillouin 
zone. (a),(c),(e),(g) Modes with fixed boundaries exhibit a 
gap. (b),(d),(f),(h) Modes with free boundaries are gapless. 
(a),(b) Despite the coupling of transverse and longitudinal 
modes, we find predominantly longitudinal and transverse 
modes on lines which we label LA (dashed blue line) and 
TA (dash-dotted red line), respectively. (c),(d) Independent 
of the boundaries, out-of-plane modes disperse quadratically 
for large wave numbers (dashed green line), (h) Free out-of- 
plane modes feature a branch with linear dispersion at the 
zone center (dashed orange line). 



a plane wave ansatz along the y direction with periodic 
boundaries. Due to their decoupling, in-plane modes 
Ux/yix,y,t) = fxfy{x)exp[i{qy — ut)] and out-of-plane 
modes Uz(x,y,t) = fz{x)exp[i{qy — cut)] can be treated 
separately. '^^ 

Exploiting the plane wave ansatz and denoting the i- 
th derivative of / as f^^\ Eq. (3) can be written as 
MxyUxJy) = -P^'^ifx, fy), where 



FIG. 4; (Color online) Displacement vector field (ux,Uy) of 
in-plane modes. Size and color of the arrows indicate the mag- 
nitude of the local deformation. We use the dimensionless co- 
ordinates X = x/W and y = y/W . (a),(b) Fixed boundaries. 
(c),(d) Free boundaries. (a),(c) Predominantly longitudinal 
modes. (b),(d) Predominantly transverse modes. 



0.2 = (l,^ g/A2), a3 = {l^iX z/g), ^ {l,i\i/q), and 

Al,2 = ±yJq^-pUJ^I{B + ll), A3,4 = ±^g2 _ ^^2/^, 

Fixed boundaries are characterized by fx{±W/2) — 
fy{±.W/2) = and by virtue of the A^, the set of linear 
equations deriving from these boundary conditions de- 
pends on the parameters q and uj. A numerical treatment 
of this linear system yields the dispersion relation [Figs. 
3;^a) and 3je)] as well as the coefficients q for the explicit 
form of the in-plane mode with fixed boundaries [Figs. 
A[a.) and 4^b)]. Other boundary conditions and the out- 
of-plane modes can be treated likewise. The eigenvalue 
problem obtained from (4) is Mzfz — (p^^/k — q'^)fzi 
where the map AAz and its eigenfunctions and eigenval- 
ues are given by 



(4) 



(2) 



Mz. fz^ f, 



and fz = Y.i=i diC^'"^ with A-^ = ±J q^ ±lu 



IV. MODE ORTHONORMALITY AND 
QUANTIZATION 



(11) 



Mxy 



-{B + fi)q^fy+^fi'^ +zBqf^'\ 



.(10) 



The general solution of this eigenvalue problem is 
{fx,fy) = X]i=i Ciajexp[Aia;], with ai = (l,zq/Ai), 



In order to quantize the vibrational spectrum of the 
graphene nanoribbon in terms of phonon creation and 
annihilation operators, the eigenfunctions of the origi- 
nal differential operators [Eqs. (3) and (4)] must be or- 
thogonal. While orthogonality of eigenmodes with differ- 
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cnt wavcnumbcrs q follows from the plane wave ansatz, 
eigenmodes with same q require orthogonal functions 
ifia,q),x, f(a,q),y) and /(a,g),^- The Indcx labels 
the phonon branch a and the wavenumber g of a spe- 
cific eigenmode. 

The map (10) is Hermitian and hence has orthog- 
onal eigenhmctions if and only if the scalar product 

!-Z/2^^{ixJy) ^xvifxJvV is real for all vector 
functions {fx,fy) in the domain of Mxy One easily 
shows via partial integration that A4xy is Hermitian if 
and only if the boundary terms satisfy 



{B + + M/;4'^ + iBqfJy _^ e R (12) 

2 

and that both fixed and free boundaries do indeed satisfy 
this condition. 

The general in-plane displacement field is 



(13) 



where the harmonic time dependence has been absorbed 
in the normal coordinate. Using the orthogonality rela- 
tions mentioned above, one can resolve the normal co- 
ordinate and derive the Lagrangian and the canonical 
momentum. The identification 



where {b(a,q)) creates (annihilates) an {a,q)- 

phonon, complies with coordinate-momentum commuta- 
tion relations, and allows for a quantum mechanical for- 
mulation of (13). Quantization of the out-of-plane modes 
is achieved in the very same way. The Hermiticity of A^^ 
follows from 



A. Elastic constants 

For graphene, most elastic constants remain to be set- 
tled by experiment and some seem to exhibit a tempera- 
ture dependence, which we do not take into account here. 
Moreover, a consistent set of constants must respect Eq. 
(9). The Zeeman energy for typical laboratory magnetic 
fields (~ 1 T) marks ~ 1 K as the temperature range 
where the phonon properties can be probed via electron 
spin relaxation. 

Cited values for Poisson's ratio a of graphene range 
from^^ 0.145 to 0.416 but accumulate around a — 0.16, 
which we use in our calculations.^^ Young's modu- 
lus of a quasi- two-dimensional material, £ = 'E^rih, fol- 
lows from its corresponding 3D bulk value and its as- 
sociated thickness h. While the most common litera- 
ture value^^"^* of EsD for graphene is ITPa, a much 
smaller value, 0.5 TPa, has been found in at least one 
experiment. We use £ = 3.4TPaA, the product of 
1 TPa and the interlayer spacing of graphite, 3.4 A. Sub- 
stituting our choices for a and £ into Eq. (9), we find 
B = 12.6eV/A2 and /x = 9.1eV/A2 for the bulk and 
shear moduli, respectively, in agreement with literature 
values."^* ^° All these values are in agreement with results 
of simulations for zero temperature.*^ 

The bending rigidity of graphene, k, is mainly deter- 
mined by the out-of-plane orbitals such that it can- 
not be inferred from other elastic constants. It has been 
shown that k decreases with increasing temperature.'*^ 
Literature values for zero temperature^^ range 
from 0.85 to 1.22 eV and we choose k = 1.1 eV. 

The mass density of graphene, p = 7.61 x lO^^kg/m^, 
follows directly from the atomic weight of natural carbon, 
12.01 u, and the interatomic distance in graphene, 1.42 A. 



B. In-plane phonons 



f* f (3) _ f *(1) f(2) _ 2„2 f* 
JzJz Jz Jz y JzJz 



(15) 



and, as above, fixed as well as free boundaries do satisfy 
this condition. The general out-of-plane displacement is 
given by 



u± = J2r{a,q)f{a,q),zeze 

a,q 



iqy 



(16) 



V. DISCUSSION OF PHONON SPECTRA 

As specific values for sound velocities, etc., depend on 
the clastic constants, wc shall first discuss these constants 
before turning to the properties of acoustic phonons. Due 
to their decoupling, in-plane and out-of-plane phonons 
can be treated separately. For each case we will consider 
fixed and free boundaries. 



The dispersion relation of in-plane modes with fixed 
boundaries is gapped and features infinitely many 
branches with different energies originating from the zone 
center. Figs. 3;^a) and 3[e). The gap relates to the 
energy necessary for fixing the boundaries and is given 
by 2.1 hy/'E/p/W. For W = 30 nm, this gap wih be 
1.0 meV, corresponding to a magnetic field of 8.4 T. For 
large wave numbers, all branches converge to a common 
line, which we label TA. A second line, labeled LA, is sup- 
ported by different branches throughout the dispersion 
relation. Due to coupling at the ribbon boundaries there 
are no purely transverse or longitudinal modes. How- 
ever, we do find that the modes on the TA (LA) line 
have predominantly transverse (longitudinal) character. 
Figs. 4(a) and 4;^b). The corresponding sound veloc- 
ities are wla = 22km/s and uta = 14km/s, indepen- 
dent of the ribbon width. These values and the ratio 
v-lk/vta = 1.6 are in good agreement with previous cal- 
culations for bulk graphene*'' (19.5 km/s, 12.2 km/s) and 
carbon nanotubes^^ (19.9km/s, 12.3km/s). Neverthe- 
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less, wc point out that our sound velocities are propor- 
tional to \pLj~p, a value that is still under discussion for 
graphene. The approach to linear, bulk-like behavior is 
expected for large wave number, where the finite ribbon 
width appears like bulk for short- wavelength phonons. 

For free boundaries, the dispersion relation of in-plane 
modes is ungapped and the two branches that start at 
zero energy converge slightly below the TA line. Figs. 
3(b) and 3(f). The sound velocities and linear behav- 
ior for large wave number do not depend on boundary 
conditions, as one would expect from the same argument 
as above. Predominantly transverse and predominantly 
longitudinal modes are shown in Figs. 4(c) and 4(d). The 
typical zero-point motion amplitude of in-plane modes is 
40 fm. 



C. Out-of-plane phonons 

The dispersion relation of out-of-plane modes with 
fixed boundaries is shown in Figs. 3(c) and 3(g). The 
gap du e to the fixed boundary conditions is given by 
22 ATi^/kTp /W^ , which yields 7.9 /ieV for W = 30 nm. 
The corresponding magnetic field is 68 mT. There are in- 
finitely many branches that correspond to different trans- 
verse excitations. Figs. 1(a) and 1(b). Again, away from 
the zone center, all branches approach bulk behavior, 
that is, a quadratic dispersion for out-of-plane modes.** 

Similarly, the out-of-plane modes with free boundaries 
disperse quadratically as in the bulk, for large wave num- 
bers. Figs. 3(d) and 3(h). The dispersion relation is gap- 
less and one branch exhibits a finite sound velocity at 
the zone center. This sound velocity amounts to about 
70 m/s for W = 20 nm, is proportional to \/k/ p /W , and 
hence goes to zero for large W, again in agreement with 
bulk graphene. The typical zero-point motion amplitude 
of out-of-plane modes is 0.4 pm. 



VI. DEFORMATION POTENTIAL AND SPIN 
RELAXATION 

Several mechanisms contribute to spin relaxation: out- 
of-plane modes via direct spin-phonon coupling and in- 
plane phonons via the deformation potential and bond- 
length changt^^. Due to inversion symmetry, piezoelec- 
tric coupling does not occur in graphene. Here, we dis- 
cuss the deformation potential, which gives the dominant 
contribution to Tj~^ under certain conditions. 

We find that any given in-plane phonon branch a cou- 
ples cither via bond-length change or via the deformation 
potential, depending on whether its displacement field is 
even or odd in the x coordinate. The branch that origi- 
nates from ujxy = 3.2 in Fig. 3 (c), labeled ag, has a flat 
dispersion at the zone center and couples to the spin only 
via the deformation potential. As a consequence of the 
density of states in Eq. ( 1 ) , this mechanism will give the 
dominant contribution if the magnetic field is tuned to a 



value where the Zeeman energy is close to the Van Hove 
singularity of ao and coupling to out-of-plane modes is 
weak. Van Hove singularities also occur for out-of-plane 
modes at different values of w^, Figs. 3(c) and 3(g). How- 
ever, u)xy and COz scale differently with W, which allows 
us to choose a ribbon width where there is a singularity 
for in-plane modes (ao) but not for out-of-plane modes. 
This situation will be discussed below. 

For the branch labeled ai in Fig. 3 (f), which is linear 
near the zone center, the spin couples to phonons only 
via the deformation potential, as well. Even though its 
density of states is finite, we discuss its contribution to 
Eq. ( 1 ) as it is in accordance with previous results for 
semiconductor quantum dots.^* 

In leading order, the deformation potential depends 
only on in-plane phonons. 



^^EPC = .goV • u^\{x,y) , 



(17) 



where go ~ 30 eV is the coupling strength^^^^ and 
V = {dx,dy). The deformation potential is independent 
of the electron spin (t,i) but it does couple different 
electron orbits {k). As a consequence, H^pc couples to 
spin indirectly when Rashba-type spin-orbit interaction, 
Hso, is taken into account. In lowest order, the spin- 
orbit-perturbed electronic states are given by 

i*t) = Mr+i;iw^»^. (18) 

where the superscript (0) indicates unperturbed product 
states. Using these spin-orbit admixed states, we find 



(fci|i/EPc|fct) 

(^EPc)fcfc'(-^So)ttfc 



(19) 



k'^k 



{HEPc)k'k{Hso)il 



+ 



Ek — + qhbB Ek — Ek' — g^isB 



where we denote the numerator in Eq. (18) as {Hso)iiJk 
and the spin-conserving transitions of i?EPC accordingly. 
This is the matrix element required to calculate the re- 
laxation rate in Eq. (1). 

We find that for a given k' the two terms in Eq. ( 19 ) ex- 
actly cancel each other at B = 0. This effect is known as 
Van Vleck cancellation and is expected for time-reversal- 
symmetric systems. Moreover, (i?so)fca vanishes if both 
k and k' are even or odd at the same time. 

For fixed GNR edges, the phonon spectrum is gapped. 
In the range 3.2 < iVxy < 3.3, the branch ao shows an 
almost fiat dispersion. Its sound velocity increases as 
0^ 1^ such that the corresponding density of states 
behaves as (pstatos)ao ^ Q^^- The matrix clement (19) 
varies as qB°-^: the dipole approximation gives rise to 
one order in q and Van Vleck cancellation* *^ *^ to one 
order in B, reduced by o;""-^ cx B~^-^ due to the prcf- 
actor in Eq. (14). In total, the contribution to the spin 
relaxation rate (1) is proportional to the magnetic field. 
Due to the Van Hove singularity of ao at ujxy = 3.2, 
we expect that T^^ (x B is the dominant behavior in 
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the range 3.8 < B < 4.0 T (12.8 < B < 13.2 T) for 
W = 100 nm {W = 30 nm), where the density of states 
of out-of-plane modes is relatively small. These are ac- 
cessible laboratory magnetic fields and hence allow for 
experimental examination of our results. 

If the magnetic field is tuned to a value where the 
Zeeman energy lies within the gap of both in-plane and 
out-of-plane phonons [Figs. 3^c) and 3^g)], the electron 
spin cannot flip due to phonon emission. Then, multiple- 
phonon processes, where the Zeeman energy corresponds 
to the difference between an absorbed and an emitted 
phonon, become important. Again due to the gap, these 
processes can be frozen out if the temperature T is low 
enough. As discussed in Sec. V, the very soft out-of- 
plane modes have a much smaller gap, which therefore 
imposes a tighter condition and which scales as W~^. 
Assuming W = 30 nm, the spin lifetime inferred from 
Eq. (1) diverges for B < 68mT and T < 90 mK. Very 
narrow GNRs with W = 10 nm are studied experimen- 
tally, as well.^^ Accordingly, the requirements for such a 
ribbon would he B < 0.61 T and T < 0.8 K. 

GNRs with free edges have ungapped phonon spectra. 
Due to energy conservation, the Zeeman energy must 
match the phonon energy, such that only the two low- 
est branches in Fig. 3 (f) are accessible for low mag- 
netic fields {B < 100 mT). The branch ai couples only 
via the deformation potential and the other branch is a 
pure shear mode. Due to its linear dispersion, we find 
B (X uj (X q and a constant density of states for ai. The 
matrix element in Eq. (19) scales as B'^'^: one order in 
B arising from each of the Van Vleck cancellation, dipole 
approximation, and the gradient in Eq. (17),^'' again 
reduced by the prefactor a B~°-^ in Eq. (14). 

Consequently, for low magnetic fields, the contribution 
of deformation potential and spin-orbit coupling to the 
spin relaxation rate (1) scales with B^ . In semiconduc- 
tors, Tf ^ DC B^ holds, as well.^^ 

VII. CONCLUSION 

Acoustic phonons are relevant for many GNR applica- 
tions and can be probed with established techniques.^" 



Using a continuum model that accounts for the 
monatomic thickness of graphene, we derive boundary 
conditions that lead to Hermitian wave equations. We 
focus on two types of boundary configurations: fixed 
and free boundaries. We explicitly give the correspond- 
ing classical solutions and, ensuring Hermiticity, infer a 
quantum theory with ribbon phonon creation and an- 
nihilation operators. Free boundaries lead to ungapped 
dispersion relations. In contrast, fixed boundaries lead 
to a gapped phonon dispersion of both in-plane and out- 
of-plane modes, which is most suitable for achieving high 
mobilities as well as long spin lifetimes. Regardless of 
the boundary configuration, all dispersion relations ap- 
proach bulk behavior for wavelengths small compared to 
the ribbon width. Sound velocities that relate to trans- 
verse and longitudinal acoustical in-plane ribbon modes 
are in good accordance with values for bulk graphene. 
We also study phonon- induced spin relaxation in GNRs. 
We find that, if the Zeeman energy is tuned close to a 
Van Hove singularity of the density of states of in-plane 
phonons, the deformation potential can be the dominant 
effect for spin relaxation. In this case, it should be possi- 
ble to probe our predicted behavior for Ti experimentally. 
If the Zeeman energy lies within the gap of both in-plane 
and out-of-plane phonons with fixed boundaries and for 
low enough temperatures, coupling to the lattice is in- 
hibited such that the spin lifetime obtained form Eq. ( 1 ) 
diverges. 
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